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ABSTRACT 

The formation of brown dwarfs (BDs) poses a key challenge to star formation theory. The observed 
dearth of nearby (< 5 AU) brown dwarf companions to solar-mass stars, known as the brown dwarf 
desert, as well as the tendency for low-mass binary systems to be more tightly-bound than stellar 
binaries, have been cited as evidence for distinct formation mechanisms for brown dwarfs and stars. In 
this paper, we explore the implications of the minimal hypothesis that brown dwarfs in binary systems 
originate via the same fundamental fragmentation mechanism as stars, within isolated, turbulent giant 
molecular cloud cores. We demonstrate analytically that the scaling of specific angular momentum 
with turbulent core mass naturally gives rise to the brown dwarf desert, as well as wide brown-dwarf 
I/"") ' binary systems. Further, we demonstrate analytically that the turbulent core fragmentation model 

also naturally predicts that very low- mass (VLM) binary and BD/BD systems are more tightly- 
bound than stellar systems. In addition, in order to capture the stochastic variation intrinsic to 
turbulence, we generate 10 4 model turbulent cores with synthetic turbulent velocity fields to show that 
the turbulent fragmentation model accommodates a small fraction of binary brown dwarfs with wide 
separations, similar to observations. Consequently, because the turbulent core fragmentation model 
can accommodate the brown dwarf desert as well as hard low-mass binaries, we suggest that these 
binary properties are not strong evidence for separate formation mechanisms for brown dwarfs. Indeed, 
the picture which emerges from the turbulent fragmentation model is that a single fragmentation 
mechanism largely shapes both stellar and brown dwarf binary distributions during formation. 
Keywords: binaries: general - ISM: clouds - stars: brown dwarfs - stars: formation - stars: low-mass 
- Turbulence 
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, 1. INTRODUCTION 

A key unresolved topic in star formation is the nature of the brown dwarf desert. The deser t refers to the observed 
relati ve absence of brown dwarf companions to solar mass stars within a distance of 5 AU ()Grether fc Lineweaverl 
120061 ) . Alternatively, Marcy and Butler 2000 described the desert by noting that 0.5% of primary stars with masses 
> 0.5Mq had brown dwarf companions within orbits of 3 AU. S everal theories have been proposed to explain the desert, 
, including both ej ection from low-order multiple stellar systems (iReipurth fe C larke 200 l|), as we ll as disk fragmentation 
( Rice et aLll2003r) followed by inward d isk migration and destruction ()Armitage fc Bonn cll 2002) or subsequent ejection 

Furthermore, it has been argued that the observed differences in the orbital properties of very low mass (VLM) 
and brown dwarf systems from those of stellar systems supports different formation mechanisms for stars and brown 
dwarfs. In addition to the brown dwarf desert, observations o f VLM and brown dwa rf binary systems reveal that these 
systems typically have narrower separations of a < 20 AU (jBurgasser et al.ll2007t ). Moreover, b inary brown dwarfs 
?H ■ typically have a significantly hi gher cutoff on their minimum binding energies than stellar systems ()Kroupa et al.l l2003. 
. . .' 120111 : iParker fc Goodwinll2011l ). In this viewpoint, a different mechanism for the formation of brown dwarfs and stars 
(e.g., disk fragmenta tion versus molecular cloud core fragmentation) may give rise to a break in the I MF across the 
subst ellar boundary (jThies fc Kr oupa 2003) as well as distinct set of binary properties for BDs and stars (|Kroupa et al.1 
120111 ). An N-body study following the evolution of very low mass binary systems, with total system mass < O.2M 
has, however, demonstrated that such tight, hard systems cannot be disrupted even in the densest stellar clusters, 
which suggests that these systems must be formed at birth ()Parker fc Goodwinll2011[ ). 

Moreover, more recent theoretical work has emp hasized limitat ions in the idealized thermodynamic treatment of 
protostellar disks in earlier star formation models (Rafikovl 120051 ). Additional theoretical models of the collapse of 
turbulent cores into protostellar disks further showed that the protostellar disks would likely remain stable to brown 
dwarf formatio n if the effects of stellar irradiation, whic h were neglected in early computational models, were properly 
accounted for ([Matzner fc Levinll2005l : lKrumholzl l2006). Subsequent radiation hydrodynamic simulations, taking into 
account stellar irradiation feedback, have shed light on brown dwarf formation process and have found that irradiated 
protostellar disks are more gravitationally stable and produce fewer brown dw arfs than models which neglected stellar 
irradiation (|Krumholz et alj|2007t iBatd 120091: lOffner et all 120091: lBatil2012D . consistent with theoretical predictions. 



State-of-the-art numerical simul ations inclu ding the effects of radiative feedback have begun to yield samples of up 
to 200 stars and brown dwarfs (Bate 2011), but cannot yet provide enough systems to make statistically-significant 
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statements about the nature of the brown dwarf desert. Consequently, the turbulent core formation model, in which 
brown dwarfs form from low -mass, turbulent-pressure confined cores, remains a viable though largely unexplored model 
(jPadoan fc Nordlundll200l . 

On the observational side, nu merous observations are consistent with the view that brown dwarfs and stars share a 
common formation mechanism ([Luhman et al.|[2007T) ■ In particular, like classical T-Tauri systems, you ng brown dwarfs 
exhib it accretion infall signatures in H-alpha and near-infrared excesses indicative of accretion disks (jMuzerolle et all 
2003). The accretion rate within these substellar systems is significantly lowe r than observed in cl assical T-Tauri 
systems, with the mass accretion rate scaling as the square of the system mass (|Mohantv et al. 1 12005). Additionally, 
brown dwarf kinematic velocities cannot be differentiated from other stars, which is consisten t with a common form ation 
mechanism for both stars and brown dwarfs and inconsistent with the ejection hypothesis (jWhite fc B asri 2003|). 

In addition, observations have highlighted the existence of a small number of widely separated, soft VLM and 
BD-BD binaries in star-forming regio ns. These syste ms include the wide young binary BD system 2MASS J11011926- 
7732 383 AB at a separation of 242 AU (lLuhmanl]2004D . wide young BD system Oph 162225-240515 at a separation of 243 



AU dJavawardhana fc Ivano v 2006; Close et al 



20071 ). wide young BD binary UScoCTIO 108 at a p rojected separation 



of 670 AU (jBeiar et all 120081 ) . and wide young BD binary FU Tau A and B (|Luh man et aT1l2009( ) at a separation of 
800 AU. In the field, obs ervations have located the VLM binary DENIS- J055146.0-443412.2 at a separation of over 200 
AU jBilleres et al.ll2005T) field VLM binary LEHPM 494/DENIS-P J0021.0- 4244 at a separation of 180 AU (|Caballerol 
2007), and VLM binary 2M0126AB at a projected separation of 5100 AU (|Artigau et al.ll2007L l200l ). 

Such widely-separated VLM and BD-BD systems are very rare and account for no more than 1% - 2% of all systems, 
but nonetheless play a crucial role in constraining models of star formation. In particular, wide separations of young 
BD/BD and BD/stellar binaries significantly wider than a typical protostella r disk size of 100 AU pose a sign ificant 
challenge for disk fragmentation models of brown dwarf binary formation (|Whitworth fc StamatelloU [2006') . The 
current VLM record-holder is 2MASS J12583501+4013083 and 2MASS J12583798+ 4014017, a very wide binary with 
both components near the BD mass limit, with a projected separation of 6700 AU (Radigan et al.l 120091 ) . Such very 
soft systems will be softened and ul timately disrupted by dynamical interactions with field stars and giant molecular 
clouds, but aslRadigan et ail (|2009D note, their estimated lifetimes, based upon standard binary evolutionary models 
(|Weinberg et al.lll987l) can still be ~ several Gyr. 

In light of the current state-of-the-art of computational simulations, a theoretical framework of a simplified model will 
help to illuminate the key physics of the brown dwarf desert and the scaling behavior of other low mass stellar systems. 
Specifically, in this pap er, we build upon and extend a turbulent fragmentation model of binary formation within 
isolated turbulent cores ([Fisher 2004) (hereafter F04), and explore the implications of this model for the formation of 
low-mass and brown dwarf binary systems. In the turbulent core model, the mass and angular momentum of binaries 
are initially established by the mass and angular momentum of an isolated, parent turbulent core. Consequently, 
the mass and angular momentum of the turbulent core gives rise to the initial, or "primordial" distributions of both 
stellar/BD and BD/BD binaries immediately subsequent to fragmentation, and prior to any subsequent dynamical 
evolution, including possible effects of competitive accretion and gravitational torques from the surrounding cluster 
protostellar members. Significantly, as we will re-derive, the turbulent core mod el predicts that the specific an gular 
momentum of turbulent cores scales with the mass of the core to the 3/4 power (Burkcrt fc Bodenheimerl l2000'). We 
will demonstrate that the scaling of specific angular momentum with turbulent core mass in turn naturally produces a 
brown dwarf desert, as well as wide brown- dwarf binary systems. Finally, we will also demonstrate that the turbulent 
fragmentation model also naturally predicts that low-mass binary and BD/stellar systems are more tightly-bound than 
stellar systems. 

We should note that the isolated core model is an idealization of the complex dynamics of giant molecular clouds. 
Historically, the concept that GMCs are highly inhomogeneous, clumpy structures with smaller, denser cores within 
these clumps was motivated by noting that if the densities of CO clumps refle cted the mean de nsity of the CMC, 
the GMC extinction would exceed that observed by over an order of magnitude (Blit z fc Shulll980l ). However, clumps 
within GMCs are far more complex than the ballistically-moving isolated objects originally envisioned by Blitz and 
Shu; current models suggest t hat they are the densest gaseous regions with in highly-complex flows governed by 
supersonic turbulent dynamics (iPadoan et al.l fl997t iPadoan fc Nord lund 20021). Moreover, recent three-dimensional 
hydrodynamical numerical simulations of supersonic turbulent giant molecular clouds demonstrate that GMC cores 
are dynamical entities, and continue to accrete from a network of parent filaments, even as their em bedded stars are 
in the process of formation (Ballcst eros-Paredes et al.ll2003t lOffner et"aT1l2010HHansen et al. 11201 lh . 

The isolated core idealization, however, remains a fruitful one because it allows us to explore the physics of tur- 
bulent fragmentation for statistically-significant numbers of binary systems. Moreover, the isolated turbulent core 
fragmentation model allows us to capture the role which tu rbulence plays in giving rise to the angular momentum of 
the parent gas distribution (Burkcrt & Bodc nheimerl l2000h . In F04, we demonstrated that this angular momentum 
distribution of parent GMC cores, in turn, plays an important role in establishing binary periods and separations. In 
particular , the isolated turbulent core model robustly agrees wi th a number of observed properties, both of prestel- 
lar cores (jBurkert fc Bodenheimerl l2000; Mat zner fc LevirJ 120051 ). and of pre-main sequence and field binaries (F04). 
Specifically, it agrees with observational data on the width of the binary period distribution. Moreover, both the model 
and observations exhibit an anticorrelation of binary period and mass ratio, and a positive correlation of binary period 
and eccentricity. The turbulent core model also predicts that low-mass model binaries originate within lower-mass 
turbulent cores with less angular momentu m, and also naturally yield narrower binaries than stellar-mass systems, in 
accord with observation (|Close et al.ll2003D . This broad agreement with the predictions of the isolated turbulent core 
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model for stellar systems motivates an extension of the model into the substellar regime. 

This paper is organized as follows. Section 2 discusses estimated scaling relationships that provide simplified predic- 
tions of our model's results. Section 3 focuses upon the methodology of our semi-analytic model for the formation of 
both stellar and substellar binaries from turbulent isothermal GMC cores, with a particular emphasis on advancements 
made to our F04 model to accommodate substellar systems. Section 4 presents the results of our model. Lastly, section 
5 presents our discussion of the results and conclusions. 

2. SCALING RELATIONS FOR TURBULENT FRAGMENTATION 
2.1. Derivation 



Turbulence in the core provides it with a net angular momentum (p3urkert fc Bodenhe imer 2000). Intuitively, as 



Landau originally suggested long ago (jLandau fe Lifshitz 1119591 iDavidsonf 2009( l. we can understand t he net angular 



momentum generated by a turbulent power spectrum consistent with Larson's linewidth-size relation (|Larsonlll981| ) 
by realizing that although the net angular momentum induced by the numerous small-scale turbulent modes will tend 
to cancel out, the few large-scale turbulent modes which fit within the core will tend to result in a small net angular 
momentum. Turbulence is inherently a stochastic physical process, and different realizations of a turbulent velocity 
field will endow a core with differing levels of angular momentum. However, the resulting net mean angular momenta 
of models of turbulent cores are n onetheless consistent with the low mean rotation ra tes implied by the linewidth 
gradients of cores mapped in NH 3 (jGoodman et al.lfl993l : iJiiina. Mvers. fe Adams Ill999h . 

Burkert and Bodenheimer 2000 showed that because most of the angular momentum endowed by turbulence is 
generated at the scale of the core, the mean specific angular moment of a population of turbulent cores could be 
reasonably estimated even under the simplifying assumption of uniform core rotation. We therefore assume uniform 
rotation to derive estimates for the specific angular momenta of model cores. Expanding upon this, we derive esti- 
mates for binary properties, including the scaling of semimajor axes, periods, and binding energies with system mass 
from the isolated turbulent core fragmentation model. These estimates capture the essential physical description of 
turbulent fragmentation, which we will subsequently elaborate upon in 331 m more detailed calculations, taking into 
account a fuller description of the stochastic variation inherent in different realizations of turbulence, as well as the 

inhomogeneous, turbulent GMC background within which individual cores are embedded. 

In our deri vation of the scaling relations, we assume that the core is a critical Bonnor-Ebert sphere (jEbert I fl 9 5 51 : 
lBonnorlll957f) with a density dist ribution modeled using an analytic lowered power law approximation 
(jNatarajan fe Lvnden-BelTI 119971 ). We calculate the moment of inertia I and the gravitational potential energy fl 
of this core as I — cM cole R 2 ole and ft = — dGM 2 orc / R corc , where c s» 0.34 and d 0.55 are constants numerically 
determined for a critical Bonnor-Ebert sphere, in terms of the core mass M coro and radius R CO rc- We take the parameter 
(5 = cR^ oic uj 2 / (2rfGM corc ) to describe the ratio of rotational energy to gravitational binding energy of the core, and find 
the specific angular momentum of the critical Bonnor-Ebert COrG (/core /-^core SC£l16S clS tllC square root of the mass and 
radius of the core; J CO re/M cmc — \ / 2cdf3GM CO r C Rcore- We combine the Larson turbulent velocity dispersion relation 
with an exponent of 1/2, a — 1.10 km s _1 .L(pc) , with the condition that the core is in virial balance; in terms of 

the virial parameter, a = 5a 2 R/(GM) a ~ 1. Consequently, we find that J col - c /M col - c cx Mclre oc Rcl? e (|Larsonlll981l ; 
ILeung. Kutner. fc Mead1ll982t iMvers I IT983). The latter scaling reflects the increase of line width with increasing size, 
and is the same scaling reported by Burkert and Bodenheimer 2000. These lead to a scaling estimate of the specific 
angular momentum with the mass of the core: 

^ = 2.6x 10- f^) V4 f^V /4 cm 2 fl (1) 

M core 1 1.3/ V0.027 V Me J 

Crucially, lower-mass turbulent cores in virial balance naturally have a lower specific angular momentum than 
more massive cores. This scaling of the specific angular momentum with mass has profound consequences for binary 
properties. While most of the mass and angular momentum in a core is carried away during the star formation process 
and does not end up in the final binary, we can describe the fractions of mass and angular momentum transferred 
from the core to the binary system in terms of a star formation efficiency, e* = M/M corc , and an angular momentum 
efficiency, ej = J/J core , where M and J are the total mass and angular momentum of the binary system. Studies 
have suggested that the star formatio n efficiency is fairly constant over a wide range of formation conditions, with a 
typical value of 0.3 (jAlves e t al. 2007). Less information is known about the angular momentum efficiency, but F04 
demonstrated that the stellar period distribution could be reproduced by a constant value of ej for a given model star 
formation efficiency e*. Thus, we use the two efficiencies to derive the scaling of the system specific angular momentum 
with total mass from the core scaling: 

M Vo.016/ V e* J VI. 3/ V0.02; \M Q J v ' 

We may then use this scaling to estimate the typical periods P and semimajor axes a of binary systems : 

'0.30\ 21/4 / a \3/4 / p \ 3 / 2 / M \ 1/4 1 (1 + q) 6 



P = 159 



0.0167 V e* J VO02/ [m^J (l - e 2)l q * d&yS [d) 
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a = 0.57(^)V^y /2 m 1/2 m (^-) 1/2 , 1 J 1+ 2 g) V (4) 
V0.016/ V e* / ^1-3/ VO-02 J \M Q J (1 - e 2 ) g 2 w 

where e is the eccentricity of the system, M 2 < Mi is the companion mass, and q = M2/M1 is the mass ratio. 

As seen above, the period and semimajor axis scale weakly with mass, to the 1/4 and 1/2 powers respectively. 
Moreover, the mass ratio of the system, q, is crucial for shaping the primordial distributions, as most of the angular 
momentum that is transferred from the core to the binary system will be associated with the orbit of the companion. 

To illustrate the significance of the mass ratio concretely, consider two 1 Mq systems : one a stellar binary with 
q = 1 and the second a BD/stellar binary with q = 0.04. For simplicity, we assume both have the mean eccentricity 
of a thermal distribution, e = 2/3. With our fiducial scalings, the semimajor axes of the stellar and BD/stellar binary 
will be 16 AU and 750 AU, respectively. The wider separation of the BD/stellar binary systems is typical of the brown 
dwarf desert. Thus, the desert naturally arises in the turbulent fragmentation model primarily as a mass ratio effect. 

Furthermore, consider a 0.16 Mq binary brown dwarf system with q — 1, which will have a semimajor axis of 8.2 
AU with our fiducial scalings. Consequently, the turbulent fragmentation model also naturally predicts that VLM 
and binary brown dwarf systems will be narrower than stellar binary systems, as they are formed within low-mass 
turbulent cores with lower specific angular momenta than stellar-mass turbulent cores. Moreover, these lower-mass 
systems will have nearer-equal mass ratios than stellar systems; this trend toward equal mass ratios further favors 
narrow-binary brown dwarf systems. 

We may expand upon our semimajor axis estimates to construct the turbulent fragmentation model prediction 
for the scaling of the minimum binding energy with the system mass. A system of mass M has binding energy of 
-Ebind = GMiM 2 /a oc M 3 / 2 (g 3 / (1 + q) 6 ). Such a binary is most weakly bound when it has a brown dwarf companion; 
these systems establish the minimum binding energies of a binary with total mass M . For a companion dwarf with 
a much lower mass than the primary, M 2 << M%, this scaling will approximately follow i5bind,min oc M -3 / 2 . As an 
example, let us consider a 0.10 Mq and a 1.00 Mq system, each with a 0.01M Q companion. We find that the former 
system is bound s» 24 times more tightly than the latter. Our result is similar to the conclusion derived by Close et al. 
2003 that VLM and binary brown dwarf systems tend to be 10 - 20 times more tightly bound that solar mass systems. 



3. METHODOLOGY 

The turbulent fragmentation model describes the formation of binary systems from turbulent giant molecular cloud 
cores. The mass and angular momentum due to turbulence in these parent cores set the mass and angular momentum 
in the resulting system, determining its orbital properties. This paper's key addition to the methodology developed in 
F04 is the relaxation of the assumption of a fixed core edge pressure. Instead, we consider cores embedded within an 
isothermal, supersonically turbulent giant molecular cloud, with a wide range of edge pressures that naturally result 
in a wide range of critical Bonnor-Ebert masses. In turn, this produces both stellar and substellar binary systems. 

Conceptually, our model of turbulent cores is similar to one recently pr oposed by (|Padoan fc Nordlundl l201lt) . 
which is based upon an earlier model of the IMF (|Padoan fc Nordlundl 120021 ). Both models assume that background 
supersonic isothermal tur bulence establishes core edg e pressures and that the core masses are uncorrelated with these 
edge pressures. However. iPadoan fc Nordlundl (|2011h consider the core mass function (CMF), as observed in a star- 
forming region, and assume that the formation of cores is distributed uniformly in time. The ongoing core accretion 
within their model leads to a local star-formation efficiency, defined in terms of the instantaneous core mass and the 
final stellar mass, greater than unity. In contrast, our model considers only the final core masses when determining 
binary properties, which are fixed by the global star-formation efficiency e». Furthermore, both models predict a 
prevalence of pressure-confined cores at substellar core masses less than the local Bonnor-Ebert mass. These cores are 
gravitationally-stable and will not form stars or brown dwarfs; thus we reject any pressure-confined cores generated in 
our model. 

Som e of the core mass and angular momentum will be lo st, through either stellar winds and outflows (jMatzner fc McKeel 
2000), or magnetic braking (|Basu fc Mouschoviaslll994Tl . respectively, during the star formation process. Thus, only 
a fraction of the original core mass and angular momentum remains in the binary star system at the time of for- 
mation. These fractions are described by two parameters, the star formation and angular momentum efficiencies, 
e* = M to t/M corc and ej = Jtot/ jeais respectively. Observational stud ies suggest that the star formation efficiency is 
e* sa 0.3, independent of mass (Al ves et alJl2007t lAndre et al.ll2010| ). This is also in agreement with Machida a nd 
Matsumoto's MHD simulation results, which suggest that cores have 0.3 < e* < 0.5 (jMachida fc M atsum otdl20"TlT ). 

To determine an angular momentum in the parent core, and the binary system through the angular momentum 
efficiency, we impose a Gaussian random turbulent field upon the core. As in F04, we set this efficiency so that the 
distribution of stellar binary systems reflects observations. For = 0.3, we find tha t e.j = 0.016 in our model. 

In F04, masses were drawn according to the Kroupa initial mass function (IMF) (|Kroupa. Gilmore. fc Toul[l99"l . 
However, this IMF had a lower cutoff at 0.08 Mq and thus exclud ed the possibili ty of modeling brown dwarfs in the 
substellar range. In this paper, we adopt the Chabrier 2005 IMF (Chabrier 2005), which extends the lower cutoff of 
the MF do wn to the deut erium-burning mass limit of 0.01 Mq. The Chabrier 2005 IMF differs from the Chabrier 
2003 IMF (|Chabrierll2003l) . which overpredicts the relative abundance of brown dwarfs to stars. By shifting the peak 
of its lognormal segment to better account for the relative numbers of brown dwarfs and stars, the Chabrier 2005 IMF 
predicts the formation of one brown dwarf between the masses of 0.03 Mq and 0.08 Mq for every four stars with 
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masses < 1 Mq. When discussing the Chabrier IMF throughout this paper, we will be referring to the Chabrier 2005 
IMF. 

The Chabrier IMF is a piecewise-defin ed function, obeying a lognormal distribution for masses < 1 Mq and a power 
law for masses > 1 Mq (jChabrierl 120051 ) . The power law regime of the Cha brier IMF has an index of —2.35 ± 0.3 
(|Chabrierll2005f ). in accordance with the famous Salpeter IMF index of —2.35 (jSalpeterlll955f ). 

To draw masses from the Chabrier IMF, we must first express it as a probability distribution function of mass 
(in Mq), to. Chabrier provides his IMF, normalized in terms of volumetric quantities, as a function of the common 
logarithm of mass, log to : 



£ (log m) 




-(log to — logTO c ) 2 /(2 x a 2 



if to < TO , 
if to > m . 



(5) 



A and B are prefactors, m c gives the location of the peak of the lognormal distribution, a describes the width of 
the lognormal distribution, —x — 1 is the power law index, and m a is the location of the break in t he IMF. Here, t hese 
parameters have values A = 0.093, B = 0.041, m c = 0.20, a = 0.55, x = -1.35 ± 0.3, and m = 1 (IChabrierll2005l). 

Co nversion between £ (log to) and £ (to) is conducted by multiplying the former by l/(TOml0) (|Scalol ll986: Ch abrierl 
2003). Therefore, the probability distribution function (PDF) of the IMF may be expressed in terms of two additional 
prefactors A* and B* as follows: 



PDF (to) 



(A*/ In 10) to 1 exp [— (logTO — logTO c ) 2 /(2 x cr 2 )] if to < to , 
(B*l In 10) m~ x ~ 1 if to > to . 



(6) 



We must normalize the PDF to unity under the constraint that the IMF is continuous at m = 1. We then calculate 
the values of A* and £?* necessary to satisfy both the continuity and normalization of the PDF : 



A* = In 10 



cr In (10) erfc 



log TO C - log TO 



aV2 



1 



■ exp 



(log TO - logTO e ) Z 
2^ 



£?* = A*m x exp 



(log TOq - logTO e ) z 
2a 2 



(7) 



(8) 



Substituting, we find that A* « 0.724 and B* rj 0.323. 

We calculate the cumulative distribution function (CDF) for the Chabrier IMF and utilize the inverse transform 
sampling method to draw masses from it for our model : 



CDF f 1 - { A * a \/^1^ erfc [( lo S TO c - l0gTO)/(crv / 2)] if TO < TOo, 

[m> ~ {C+ (B*/lnlQ)(m- x /x)(l - (m/m )~ x ) if to > m . 



(9) 



Here C « 0.896, the value of the CDF at to = 1 Mq, ensures the continuity of the function across the break in masses. 

We assume that binary masses are uncorrelated and thus draw both masses in each system independently from the 
IMF. This assumption works well for lo w stellar masses, but breaks down at high masses; O and B type stars tend 
to be paired with near equal mass stars (San a et aLll2009l) . Therefore, we account for the limits of random pairing by 
restricting all stars in our model to < 2 Mq. Furthermore, by drawing the masses as an input in this manner, our 
model assumes that each of its stars and brown dwarfs primordially form as part of a binary. Additionally, as our 
model is primordial, we will not directly capture any dynamical evolutio n of the binary popula tion as cluster evolution 
causes soft binaries to become softer and hard binaries become harder ((Weinberg et al.lll987[ ). 

Once we have drawn the primary and companion masses, we set the parent core mass by M core = M tot /e*, where 
M to t is the total binary system ma ss. We assume that the core is described by a turbulently-supported Bonnor-Ebert 
sphere (|Ebert 1119551 IBonnorlll95"7l ) at or exceeding the critical edge density at which the core becomes Jeans unstable. 
Once we determine a core edge density, we utilize a Bonnor-Ebert sphere density profile to provide the core's overall 
density structure. Combining this information with a Gaussian random turbulent velocity field, we may calculate a 
realization of the core angular momentum. 

We assume a C MC model with mass M = 10 5 Mq in virial equilibrium. To do so, we take the 2500 Mq GMC 
model assumed bv lKrumholz fc McKcc (2005) and scale it to 10 5 Mq under an isothermal temperature of 10 K. Using 
this model GMC, we calculate a dispersion velocity of er c ; = 4.0 km s _1 and then use the mean density relationship 
derived by Krumholz and McKee 2005 to find an mean mass density, p c \: 
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(10) 



From this, we derive a mean model GMC mass density of 4.8 x 10 g cm . Using standard mass fractions of 
molecular hydrogen and helium, this corresponds to a number density of ~ 1265 cm -3 , which provides the expectation 
value of our lognormal distribution. This allows us to draw the core edge densities from the distribution. 



G 



Jumper, P.H. and Fisher, R. T. 



Supersonic isother mal turbulence produces a lognormal distribution of pressures (and thus densities) within GMCs 
(|Padoan e t al. 1997). Thus, we draw the edge densities of our model cores from such a distribution. However, as we 
have drawn our stars and brown dwarfs from the IMF, we have effectively assumed that the cores in our model must 
collapse. We accept edge densities that make the cores Jeans-unstable; any lower densities resulting in Jeans-stable 
cores are rejected and redrawn. 

Padoan and Nordlund give a log-normal probability distribution function for number density n: 



p(n)dn 



1 



(27TCT 2 ) 1 /2 n 



exp 



1 / In (n) - In (n) 

2 I a 



(11) 



based on an mean number density of 1 cm , which results in Inn = —o /2 (jPadoan &: Nordlund! l200l . For different 
mean densities, the relationship between the mean of the distribution [l and Inn becomes Inn = ln(/i) — a 2 /2. We 
utilize the mean number density determined previously (~ 1265 cm -3 ), substitute this into the number density PDF 
and utilize the inverse transform sampling method as before to draw edge densities from the CDF of the Padoan and 
Nordlund log-normal distribution. 

Padoan and Nordlund describe a 2 in terms of two parameters, b, which reflects the degree to which the model's 
turbulent forcing is solenoidal or compressive, an d Ms, the turbulent Mach number, such that a 2 = ln(l + b 2 M 2 .) 
(jPadoan fc Nordlund! I2004t iFederrath et aT1l2010f). Models a ssuming solendoial forcing have included proposals of 
b = 0.26 (|Kritsuk et al. II2007T ) and b = 1/3 (jPrice et al.|[20Tl[) . The actual nature of the forcing is unk nown, but it is 
likely an intermediate between solenoidal and compressive, similar to the original proposal of 6 = 0.5 (jPadoan et al.l 



3iy 

I1997T) . which we will adopt for our model. 

In addtion to both of the binary system's masses, the parent core's edge density, and Gaussian random turbulence, 
we must also draw the eccentricity of the final binary system to complete our orbital calculations. This value is drawn 
from the thermal distribution, /(e) = 2e, which has a mean value of 2/3. Together, these values yield the binary's mass, 
angular momentum, eccentricty, and mass ratio, which are then used to calculate the binary's period and semimajor 
axes: 



,'2tt 



J 3 



l 



(i + <z) e 



M 5 J (1-e 2 )! 



(12) 



a= G\M 



1 1 (1 + g) 4 

(1 - e 2 ) M q 2 



(13) 



We repeat these procedures for each binary, generating 10 4 model systems. 



4. RESULTS 
4.1. The Primordial Brown Dwarf Desert 

The key feature of the brown dwarf des ert is the relative absence of brown dwarf companion s to stars comparable to 
a solar mass within distances of 3-5 AU (jMarcv fc Butlerl 120001 : IGrether fc Lineweave r 2006). Studies have proposed 
upper limits between . 5% and 1.0% on stars exceeding 0.5 M w havin g brown dwarf companions within these distances 
(|Marcv fc Butlerll2000l: iGizis et alJl200H IGrether fc Linewe aver 20061) . In our model, we drew 3709 primary stars with 
M > 0.5M Q out of a sample of 10 4 systems. Of these systems, 7, or 0.19%, had brown dwarf companions within 3 
AU of the primary, while 12, or 0.32%, had such companions within 5 AU. Thus, turbulent fragmentation naturally 
accommodates the brown dwarf desert. 

We further demonstrate the dearth of narrow brown dwarf companions with Figure[TJ which describes the distribution 
of semimajor axes for all model brown dwarf companions to primary stars between 0.5 Mq and 2 Mq under the 
efficiencies of e* = 0.30 and ej — 0.016, a sample of 972 systems. Among this population, narrow systems remain 
scarce. Only 21 of 972 such systems, or 2.2%, have separations < 10 AU, while 166, or 17%, have separations 
< 100 AU. 

Moreover, the peak of the semimajor axis distribution is located on the order of 1000 AU, as seen in Figure [TJ 
demonstrating our model brown dwarfs' preference for wide separations. This is further supported by additional 
measures of central tendency, including a median of the semimajor axes of ss 911 AU. Thus, our results are consistent 
with the conclusions of Gizis et al. 2001, who suggested that brown dwarf companions may be common at distances 
greater than 1000 AU. 

Crucially, the current model only considers the primordial orbital parameters that result from the initial formation 
of its binary systems. Therefore, the production of a brown dwarf desert suggests that turbulent fragmentation may 
explain the desert primordially, prior to the influence of further evolutionary effects. This contrasts with the previous 
evolutionary models explaining the emergence of the brown dwarf desert from later dynamical effects, such as migration 
or ejection, as noted in the introduction. Such evolutionary mechanisms continue to sha pe the binary distribution well 
past formation, but the binary distributions themselves are largely established at birth ([Parker fc Good win 2011). 
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4.2. Brown Dwarf - Brown Dwarf Binaries 

Binary brown dwarf systems preferentially form at narrow separations. In a review paper, Burgasser et al. 2007 
reported that 93% of VLM systems have separations less than 20 AU. Additionally, they indicated that the peak of the 
binary brown dwarf semimajor axis distribution is located between 3 AU and 10 AU, with a mean of 4.6 AU. Similarly, 
Close et al. 2003 reported that the peak of this distribution was located at sa 4 AU. Our model produces similar 
results, as shown in Figure [2l which provides data for 567 binary brown dwarfs systems at e* = 0.3 and ej = 0.016. 
We find that 512, or 90.3% of our model systems have separations less than or equal to 20 AU, with 469, or 82.7% of 
systems within the narrower limit of 10 AU. Our distribution appears to peak at a somewhat lower value than those 
reported by Burgasser et al. and Close et al, peaking in the vicinity of its median of « 1.63 AU. 

Furthermore, the mean of our distribution, « 15.4 AU, initially appears to contrast with the observational results of 
Burgasser, Close, and coauthors. However, due to limited statistics, most observational samples do not capture many 
of the widest binary brown dwarf systems, which may have semimajor axes in the hundreds to thousands of AUs, with 
a current record 6700 AU (Rad igan et al.ll2009h . In contrast, our model has captured a long tail in the primordial 
semimajor axis distribution, including 12 systems with separations in excess of 100 AU, extending to a maximum of 
a R3 3481 AU. This tail results from turbulence in cloud cores; different realizations of the turbulence in the the same 
core may result in greatly different angular momenta, and in turn orbital separations. This inherent stochastic nature 
of turbulence gives rise to unusually wide binary brown dwarfs, which increase the mean separation of our sample. 
Moreover, the primordial nature of the tail is significant, as soft, or weakly bound, binaries tend to become softer 
and more widely separated over time. Therefore, while would not expect to find a primoridal separation rivaling the 
current record holder in the field, we could expect that a time evolution of the model systems would produce more 
wide systems and some systems with separations of a few thousand AU. Thus, the turbulent fragmentation model can 
accomodate both the existence of a thin tail of wide systems as well as the more typical, narrower systems. 

4.3. Binding Energies 

Observational studies have found that very low mass and binary brown dwarf systems tend to be bound more tightly 
than stellar systems, demonstrating that the minimum binding energy for VLM a nd BD binaries is roughly 10-20 times 
higher than that for solar mass systems (|Close et al.l l2003; Bur gasser et al.| [2007). In Figure [31 we compare our model 
results against a compilation of previous observational results presented in Burgasser et al. 2007's Figure 6. Circles 
represent our model systems in Figure[3j while stars represent the observational results. Each set is represented by 191 
systems. Additionally, we draw a dashed line representing the approximate scaling of the minimum binding energy 
derived in §2, -Ebind.min oc Af~ 3 / 2 . We base this line from the binding energy predicted for a 1.01 M Q system with a 
1.00 Mq primary and a 0.01 M Q secondary (q — 0.01) by our scaling estimates, 2.9 x 10 40 ergs. 

As seen in Figure[3l our model results produce an trend of increasing minimum system binding energy with decreasing 
total system mass, consistent with the expectation presented by observation. Both the observations and our model 
include a few far more weakly bound outliers. Moreover, the minimum binding energies presented by both our data 
points and our scaling estimate line are similar to those found observationally and to each other. More generally, there 
is a wide overlap of the model and observational results throughout the parameter space. Ultimately, the similarities 
betweens our model and the observational results suggest that the physics of turbulent fragmentation during primordial 
formation may play a key role in establishing the observed trend in minimum binding energies. 

5. CONCLUSIONS 

The turbulent fragmentation model provides an alternative to the hypothesis that brown dwarfs form a distinct 
population with a separate formation m echanism from stars and that the IMF is discontinuous across the substellar 
mass boundary (jThies fc Kr oupa 2007]). In particular, the brown dwarf desert, the greater binding energies of VLM 
and BD/BD systems compared to solar mass systems, and the result that "one BD is pr oduced for every 4-6 fo rmed 
stars" have all been used as evidence to argue in favor of such distinct populations (jThies fc Kroupal 120071 ) . In 
this paper, we have demonstrated that, when properly scaled to account for both the turbulent linewidth as well as 
turbulent core edge pressure, the predictions of turbulent fragmentation can accommodate observed systematic trends 
in binary properties. Therefore, given that a turbulent fragmentation model may also explain these binary brown 
dwarf properties in terms of a single core fragmentation mechanism, we suggest that these binary properties are not 
strong evidence for separate formation mechanisms for brown dwarfs. Indeed, the picture which emerges from the 
turbulent fragmentation model is that a single fragmentation mechanism largely shapes both stellar and brown dwarf 
binary distributions during formation. Well-understood, subsequent dynamical interactions, b oth within the nascen t 
stellar cluster, and with field stars and GMCs will continue to evolve the binary distributions ()Weinberg et al.|[l987f) . 
but these interactions alone cannot account for the observed binary distributions of VLM and BD binaries, as recent 
N-body models have demonstrated (P arker k. Goo dwin 20Tl(). 

One deficiency of the turbulent fragmentation model presented here is that it assumes equal binarities for stars 
and brown dwarfs, whereas observations have demonstrated that the overall field binarities for VLM stars and brown 
dwarfs are lower than those of G dwarfs (jDuquennoy fc Mavoriri991l : iFischer fc Marcvlll992t iClose et all 120021 ) . Our 
current model addresses only the primordial binary distributions and consequently lacks sufficient depth to address this 
issue in its entirety. However, a general outcome of turbulent fragmentation is that the most weakly-bound systems 
preferentially include those with brown dwarf or VLM companions. Consequently, we conjecture that the softening 
and eventual disruption of these most weakly-bound systems will result in a decreased binarity fraction from stellar 
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to brown dwarf masses. Other authors, evolving VLM systems, have com e to a similar conclusion that the brown 
dwarf binary fraction must have been higher at birth than is now observed (iParker fc Go odwin 20lT[). There is some 
evidence for this trend in Figure 3, which demonstrates that several brown dwarf binaries are significantly softer than 
turbulent fragmentation predicts for the primordial distribution. These ultrasoft BD binaries in the field may be the 
result of subsequent softening of initially soft primordial BD binaries. However, more work on the evolution of these 
binaries in larger clusters is necessary in order to quantify the evolution of the binary fractions themselves. 

One key future test of the turbulent fragmentation model predictions for brown dwarf binaries will be the development 
of numerical simulations of giant molecular clouds in virial equilibrium which can produce statistically large (N > 10 4 ) 
samples of binaries. In addition to allowing for a closer comparison to observations in the field, such models may allow 
further insight into the problems of the BD and VLM binarity fractions and the relative frequencies of BD/BD and 
BD/stellar systems. Such simulations will be challenging, but may be within the capability of the next-generation of 
petascale supercomputers. 

In addition, the turbulent fragmentation model prediction that low-mass stellar and brown dwarf binaries will have 
a higher binding energy than solar mass binaries may have broader applicatio n to our understanding of the core mass 
function and the origin of the IMF. In particular, iPadoan fc Nordlundl ([201 1[ ) have recently described how supersonic 
isothermal turbulence shapes the observable pre-stellar phase of turbulent low-mass cores in principle. In practice, 
however, it will be extremely challenging to distinguish in both observation and simulation the truly pre-stellar cores 
from the numerous unbound turbulent transient density fluctuations, which dominate the core mass spectrum at very 
low masses. However, because the lowest-mass binaries also fragment directly from the lowest-mass turbulent cores, the 
turbulent core fragmentation model suggests that the turbulent pre-stellar core phase can be traced in well-established 
properties of the stellar binary field distribution. This connection between turbulent cores and binaries may thereby 
allow a direct and relatively clean window into the physics which shapes the turbulent core mass spectrum and, in 
turn, the initial mass function. 

The authors acknowledge discussions with David Ribeiro. RF acknowledges research support from NSF Grant CNS- 
0959382 and AFOSR DURIP Grant FA9550-10-1-0354. PJ acknowledges support from NSF Grant DMS-0802974. 
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Table 1 

Distribution of parameters for all stellar systems without brown dwarfs at e* = 0.3. As in F04, ej is set to fit the median log period of 

stars. 
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Figure 1. Probability distribution of log semimajor axes for 972 model stellar/BD binary systems with primary masses > 0.5Mq, e* = 0.3, 
and ej = 0.016. Error bars have been calculated with Poisson statistics. 
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Figure 2. Probability distribution of log semimajor axes for 567 model BD/BD binary systems with e* = 0.3 and ej = 0.016. Error bars 
have been calculated with Poisson statistics. 
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Figure 3. Binding energy versus total system mass for 191 observed systems from data compiled by Burgasser et al. 2007 (stars) and a 

commensurate number of model systems (circles). The dashed line represents the approximate scaling of i?bind,min M -3 / 2 . The solid 
line represents the minimum mass model binary system. 



